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Explicit Robinson-Trautman solution with minimally coupled free scalar field is derived and 
analyzed. It is shown that this solution contains curvature singularity which is initially naked but 
later the horizon envelopes it. We use quasilocal horizon definition and prove its existence in later 
retarded times using sub- and supersolution method combined with growth estimates. We show 
that the solution is generally of algebraic type II but reduces to type D in spherical symmetry. 
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I. INTRODUCTION 

Solutions to Einstein equations with scalar field source 
provide very useful tool for understanding relativity due 
to the simplicity of the source. Recently, it becomes 
progressively plausible that such fields might really exist 
(LHC) and potentially play a fundamental role in physics. 
In classical General Relativity they were used to study 
counterexamples to black hole no-hair theorems and in 
many other areas. These results were mostly based on 
highly symmetric solutions and it is therefore important 
to provide solutions with less or no symmetries to sub¬ 
sequently analyze if those results hold in more generic 
situations and are not tied to a specific symmetry. 

Robinson-Trautman spacetimes represent an impor¬ 
tant class of expanding nontwisting and nonshearing so¬ 
lutions GH describing non-spherical generalizations of 
black holes. In general, they do not posses any Killing 
vectors thus providing important solutions devoid of sym¬ 
metry. Many properties of this family in four dimensions 
have been studied, especially in the last 25 years. In par¬ 
ticular, the existence, asymptotic behaviour and global 
structure of vacuum Robinson-Trautman spacetimes of 
type II with spherical topology were investigated by 
Chrusciel and Singleton |4|-[6j. Robinson-Trautman so¬ 
lutions were shown to exist for generic, arbitrarily strong 
smooth initial data for all positive retarded times, and to 
converge asymptotically to corresponding Schwarzschild 
metric. Extensions across the “Schwarzschild-like” fu¬ 
ture event horizon can only be made with a finite order 
of smoothness. These results were generalized in 0, H] to 
Robinson-Trautman vacuum spacetimes with cosmologi¬ 
cal constant. These cosmological solutions settle down to 
a Schwarzschild-(anti-)de Sitter solution at large times 
u. Finally, the Chrusciel-Singleton analysis was ex¬ 
tended to Robinson-Trautman spacetimes including mat¬ 
ter, namely pure radiation [9j. IQ], showing that they ap¬ 
proach the spherically symmetric Vaidya-(anti-)de Sitter 
metric. Generally, the solutions of this family settle down 
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to physically important solutions. The location of the 
horizon together with its general existence and unique¬ 
ness for the vacuum Robinson-Trautman solutions has 
been studied by Tod [11]. Later, Chow and Lun [13 
analyzed some other useful properties of this horizon 
and made numerical study of both the horizon equa¬ 
tion and Robinson-Trautman equation. These results 
were later extended to nonvanishing cosmological con¬ 
stant [H|. The anisotropy of Robinson-Trautman hori¬ 
zon and its associated asymptotic momentum was also 
used in the analytic explanation of an "antikick" appear¬ 
ing in numerical studies of binary black hole mergers [l5j . 

Robinson-Trautman spacetimes (containing aligned 
pure radiation and a cosmological constant) were also 
generalized to any dimension 11 ill. Existence of horizons 
was subsequently analyzed in [161 ]. Finally, Robinson- 
Trautman solutions with p-forrn fields in arbitrary di¬ 
mension were derived recently G3- One of the re¬ 
sults mentioned therein rules out the existence of aligned 
scalar field (where alignment refers to the gradient of the 
field) for generic Robinson-Trautman case. 

The solutions for "stringy" Robinson-Trautman space- 
time corresponding to Einstein-Maxwell-dilaton system 
were obtained in 18]. Recently, scalar field solutions for 
Einstein- Maxwell-Lambda system with a conformally 
coupled scalar field belonging to Plebanski -Demiahski 
family (containing type D solutions of Robinson- 
Trautman class) were derived in [l9j - 

II. VACUUM ROBINSON-TRAUTMAN 
METRIC AND FIELD EQUATIONS 

The general form of a vacuum Robinson-Trautman 
spacetime can be given by the following line element GJ- 

M 

r 2 

ds 2 = —2 H d u 2 — 2 dit dr + (di/ 2 + da; 2 ), (2.1) 

P 2 

where 2 H = A( In P) — 2r(InP) iU — 2 m/r — (A/3)r 2 , 

A =P\d XX +dyy), (2.2) 

and A is the cosmological constant. The metric depends 
on two functions, P(u,x,y) and m(u) , which satisfy the 
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nonlinear Robinson-Trautman equation 

AA( In P) + 12 m( In P), u — 4 ?n iU = 0 . (2-3) 

The function m(u) might be set to a constant by suitable 
coordinate transformation for vacuum solution. 

The spacetime admits a geodesic, shearfree, twistfree 
and expanding null congruence generated by k = d r . The 
coordinate r is an affine parameter along this congruence, 
u is a retarded time coordinate, and x, y are spatial coor¬ 
dinates spanning transversal 2-space with their Gaussian 
curvature (for r = 1) being given by 

IC(x, y, u) = A(lnP). (2-4) 

For general fixed values of r and u, the Gaussian cur¬ 
vature is /C/r 2 so that, as r —> oo, they become locally 
flat. 


III. SOLUTION COUPLED TO A SCALAR 
FIELD 

We consider the following action, describing a scalar 
field minimally coupled to gravity, 

s = J (Px^r^in + v^vv] (3.i) 

where 1Z is the Ricci scalar for the metric g pl/ . The mass¬ 
less scalar field ip is supposed to be real and we use units 
in which c = H = 8iiG = 1. By applying the variation 
with respect to the metric for the action m, we get 
Einstein equations 

Put' — -^g^vTZ = T (3.2) 

The energy momentum tensor generated by the scalar 
field is given by 

T a ipS7pip (3.3) 

and the scalar field must satisfy corresponding field equa¬ 
tion 


Oip(u,r) = 0 (3.4) 

where □ is a standard d’Alembert operator for our metric 

(13.51) . 

For the matter of convenience we will be looking for 
the metric in the following form 

ds 2 = —2(H(u,r) + K(u,x,y))du 2 — 2dudr 

+ p[ M ’ r L (dx 2 + dy 2 ) (3.5) 

P(x,y) 


of E3- The nontrivial components of the Ricci tensor 
corresponding to the metric (13.51) are 


1Z UU = 2(2^H, r + H, rr ^(H + K)+2^(H + K) <u 

2 P 2 

— + R,uu ) + -^2 {K,xx + R,yy) 

Krr = -2^ (3.6) 

K 

It T* H T* it 'Vll 

IZru = ^ur — 2 — 11 -b H.rr 

K 

P-xx = Pyy = — y) + 2(H + -R)(PP,r),r + 

+2 RR r H >r — 2 (RR u ) }r } 

where as usual () iX . = ^r() and 

k(x,y) = A(ln P(x,y)) (3.7) 


where A is still given by expression (12.21) with P replaced 
by P. 

We will use the following form of equations equivalent 
to Einstein equations (|3.2D coupled to energy momentum 
tensor (El 




( <P 2 u <Au¥V 0 0 \ 

0 0 

o o o o 

\ 0 0 00 / 


(3.8) 


From the above equations describing gravitational field 
and field equation for the scalar field (|3.4D we obtain the 
following expressions for unknown metric functions and 
scalar field 


H(u, r) 
R(u, r) 
K(u, x, y) 

<p(u,r) 

A k(x,y ) 
U(u) 


dU(u ) 


2 U(u) du 


lU(t 


U{u) 


Hx,y) 

2 U(u) 

V2 


U(u)r - Cp 
U(u)r + Co 


ye 


(3.9) 


in which Co ^ 0, a, 77 , 7 , u are constants and ui = In 
the following we will assume Co > 0 ,a> 0 ,? 7 > 0 , 7>0 
for simplicity of discussion. 


IV. PROPERTIES OF THE SOLUTION 


The scalar field is assumed to be function of u and r only 
(ip(u,r)). The dependence on r means that the scalar 
field is not aligned and thus is not ruled out by the results 


First, we should ensure that our solution really be¬ 
longs to the Robinson-Trautman family. This is simply 
confirmed by studying the properties of a null congruence 
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generated by vector 1 = d r . Such congruence is geodesic, 
nontwisting, nonshearing and its expansion is given by 


n R r 2 U(u) 2 r 
1 R U(u) 2 r 2 - C 2 ' 


(4.1) 


Evidently, the above expression is positive only for r > 
which may seenr not satisfactory. However by in¬ 
specting the Kretschmann scalar 


1 

R(u, r ) 8 


(4.2) 


and using (13.91) we immediately see that the geometry has 
singularities for r = ±ro = ± J)°, . Naturally, we are led 

to constrain the range of coordinates to r £ ^^^y,oo^. 
In this range the expansion (14.11) is everywhere positive, 
diverges at the singularity and approaches zero at infinity 
(as r —> oo). Also, one can check from the line element 
that the singularity is a standard pointlike one. Due to 
the asymptotic behaviour of function U(u) (see (13.91) ) the 
singularity tends to ro = 0 as u —> oo. The singularity 
appears due to the divergence of the scalar field and its 
energy momentum tensor. 

Asymptotically (u —> oo), the scalar field itself is van¬ 
ishing everywhere outside the singularity (see (13.91) ) while 
it diverges at r = 0. So there would be no scalar hair 
left outside when the spacetime settles down to the fi¬ 
nal state. Indeed, our geometry approaches the original 
Robinson-Trautman form for u —^ oo when we de¬ 
fine P(u,x,y) = P(x,y)/U(u). In this case one can ap¬ 
ply the Chrusciel-Singleton analysis of asymptotic 
behaviour to recover the spherical symmetry of the final 
state which neccesarrily points to Schwarzschild solution. 

When the singularity is present in our solution we will 
investigate if it is covered by a horizon. Due to dynamical 
nature of the spacetime it is preferable to use the quasilo¬ 
cal definitions of horizon — apparent horizon [21J , trap¬ 
ping horizon [22| or dynamical horizon [23| . The basic 
local condition is shared by all the standard horizon def¬ 
initions: these horizons are sliced by marginally trapped 
surfaces with vanishing expansion of outgoing (ingoing) 
null congruence orthogonal to the surface. We will be 
looking for the horizon hypersurface in the following form 


r = A4(u, x, y) 


(4.3) 


and study the expansion of compact slices of such hyper¬ 
surface given by u = uq = const, (with JVi(uo,x,y) = 
M(x,y)). The requirement of compactness necessarily 
means that the two-spaces spanned by x and y are com¬ 
pact as well. We construct null vector fields orthogonal 
to surface r = M(x, y) 


1 = 8 r 
k = d u 


(4.4) 


-(M 2 + M 2 )-(H + K) 


+-^2 {M x d x + Af y 9y) (4-5) 


that satisfy normalization condition 1 • k = —1. From 
the geometry of the situation one can deduce that con¬ 
gruence 1 is outgoing while k is ingoing. The expansion 
of the congruence generated by 1 is always positive (see 
m and the discussion beneath) so we are looking for 
vanishing of expansion related to the other congruence 
k. These two conditions (0i > 0 and 0i = 0) mean 
that we are looking for the past horizon according to the 
definition given by [22j. The second expansion is given 

by 


0 k = [AM — (Ini?) r (VM • VM)— 

-(K + H)(R 2 ) >r + (R 2 ) !U ], (4.6) 

where Laplace operator and scalar product denoted by 
dot correspond to metric hijdx l dx ; l = p ( X|j/ )2 ( dx 2 +dy 2 ) 
on the space E spanned by x,y. So the horizon is given 
by the solution of the following quasilinear elliptic partial 
differential equation 

{AM - (In R), r (VM ■ VM) - (K + M)(R 2 ) >r 

+ (-R 2 ),u} | r=M(x,y)ku=uo = 0 (4-7) 

where all dependence on r is replaced by the function 
M(x,y) and u is evaluated to arbitrary constant value 
Mo¬ 
lt is impossible to solve this equation generally but 
fortunately we can get some useful information about 
the existence of solution using the technique developed 
for the case of Robinson-Trautman spacetime in higher 
dimensions jl6||. The proof of existence of the solution to 
the same type of quasilinear equation (Ait = F(x , it, Vit)) 
was given there by combining several steps motivated by 
[24l | and using results from [25l - [27j . The main issues were 
to provide an estimate for the function F of the form 
\F\ < B(u)(l + |Vit| 2 ) (where B{u) is increasing function 
on K + ), and to show the existence of a sub- and a super¬ 
solution (28| u~ < u + , £ C 1,/3 (E) D L°°( E) (here 

C 1, P( E) are Holder continuous functions of some suitable 
index /3). Then we know there is a solution u £ C 2 ' L ( E) 
(for some i) satisfying u~ < u < it + . 

In our case, to provide an estimate of the form 
|M| < B{M){ 1 + |VM| 2 ) (the norm is taken with re¬ 
spect to the two-dimensional metric hij) for the hori¬ 
zon equation (1X71) when considered in the form AM = 
F(x, y, M, VM) where 


F = (In R(uq, M)) r |VM| 2 + k(x, y)M - 

U 2 (u 0 ) 

(4.8) 

one has to deal with the singular behaviour of (lni?) jr 
at r = AA \y e can d 0 this either by removing 
the vicinity of singularity from our domain r £ R + \ 
~ 5 )> j^y(l +<5)^ or by continuing (with some 
appropriate smoothing) the divergent function on the 
problematic interval f ~ 5), + S)) with a 
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constant value it attains on the boundary of the inter¬ 
val. Now, with all the coefficients of the equation finite 
one can construct the bounding function B(u) easily and 
thus we can proceed to the construction of sub- and su¬ 
persolutions M±. 

First, we note that due to the selection of sign for the 
free constants made at the end of previous section (Co > 
0, a > 0, r] > 0, 7 > 0) we obtain U jU > 0 if we restrict 
our attention to retarded time region u £ We 

can then understand our solution as being given by ini¬ 
tial conditions specified at Ui n = — which corresponds 
to usual understanding of Robinson-Trautman solution. 
As usual, we are looking for constant sub- and superso¬ 
lutions but we are unable to provide them independently 
of the value of uq. Generally, we can find the sub- and 
supersolutions in the following cases: 

... ^ ... _ min(fc(x,j/))—max(fc(x,y))<5— C 0 ri 

• u 0 < m - -- 

M~ = 0 (4.9) 

M+ = mk {1 - s) (4 ' 10) 




„. ^ „. _ max(k(x,y))(l+5)-C 0 V 

U0>U 2 - -2C^A- 


M~ 

M+ 


Co 

U(u 0 ) 


(1 + - 5 ) 


CjU u (u 0 ) 
min(fc(x, y))U 2 (u 0 ) 


( 4 . 11 ) 

( 4 . 12 ) 


V. ALGEBRAIC TYPE OF THE SOLUTION 


Now, we would like to see if the geometry of our space- 
time is sufficiently general. Since vacuum Robinson- 
Trautman spacetime is generally of algebraic type II we 
would like our solution to be at least of the same type 
and not more special. Our preferred tetrad for determin¬ 
ing the Weyl scalars of our solution is given by different 
null vectors compared to (14.41) 


1 = d r 

k = 8 U - (H + K)d r ( 5 . 1 ) 

m = ~^R ydx + IdV ' > 


where I is a complex unit. The Weyl spinor computed 
from this tetrad has only the following nonzero compo¬ 
nents 


’Fo 


’Fi 

*F 2 


1 


^ P{k,yy k,xx + 


AUR 2 

-(fc,X - /fc y )(P j2 . - IP >y 
V2PR 


UR 2 
1 


( I k,y ) 

[Uk - ( U, u r + k)(RR, rr R,r 2 y 


6 UR 2 

-2 URR,ru + {RU iU + 2 UR^ r )R tr \ 


( 5 . 2 ) 


Both bounds u\ and u 2 are in the restricted range of co¬ 
ordinate u. Evidently, the first case would provide exis¬ 
tence of solution only beneath the position of singularity 
(or, in other words, inside the singularity) which is irrele¬ 
vant and moreover we have already restricted the range of 

r £ (win y, 00 )- I n the second case, one can easily check 

that the necessary condition M~ < M + is indeed satis¬ 
fied for uq > u 2 and we certainly have a horizon given by 
r = M{x,y) where M~ < M(x,y ) < M + . Note that we 
suppose that min(/e(a:, y)) > 0 for the last estimate M + 
to be valid. 

If we allow k m i n = min (k(x,y)) < 0 (we define ac¬ 
cordingly k max = max(fc(i, y))) we are unable to provide 
constant supersolution in the case uq > u 2 . Instead we 
can use the knowledge of how Laplace operator acts on 
k(x,y) (13.91) and the observation that the first term of 
the definition of function F (TO) is always positive to 
provide the following non-constant supersolution 

M + = c[k max - k(x, y )] + (1 + 6 ) (4.13) 

U(uo) 

where 

c _ C 0 [C 0 U, u (u 0 ) - k min U(u 0 )(l + (5)] ^ ^ 

U^(v>o){cP T kmaxkmin k m in ) 

This estimate works if (or + kmaxkmin — kC n ) > 0. 

Even-though our solution possesses singularity at any 
retarded time u this singularity appears to be initially 
naked and the horizon develops only in later time. 


Now, we can easily determine the type irrespective of 
possible non-optimal choice of tetrad by using the review 
of explicit methods for determining the algebraic type 
in m that are based on 3(1 . Namely, when we use 
invariants 

/ ^f 4 «k 2 \ 

/ = >F 0^4 - 4^1 ^3 + 3 ^ 2 ! J = det I ^2 ’Fi 

\ ^2 'Fl ’Fo J 

we can immediately confirm that I 3 = 27 J 2 is satisfied so 
that we are dealing with type II or more special. At the 
same time generally IJ ^ 0 so it cannot be just type III. 
Moreover, the spinor covariant Rabcdef has nonzero 
components 

Poooooo = ’Fi(3’I , o'I , 2 — 2'F 4 ) (5-3) 

-R 000001 = - 4' 2 (3'I' o 'I ' 2 — 24 / 3 ) (5-4) 

which means that generally the spacetime cannot be of 
type D. So indeed our scalar field solution is of the most 
general type possible for the Robinson-Trautman vacuum 
class. Which does not mean that there cannot be a scalar 
field solution of type I. Moreover, inspecting the com¬ 
ponents of the Weyl spinor (15.21) one concludes that in 
the special case of k{x,y) = const > 0 (constant posi¬ 
tive Gaussian curvature of compact two-space spanned by 
x, y) the algebraic type becomes D consistent with spher¬ 
ical symmetry. Finally, since U/i = 0 implies Tq = 0 we 
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cannot have all components of spinor covariant Qabcd 
(see [Hand [H) vanishing while having nonvanishing 
Weyl spinor. This means that our family of solutions 
does not contain type N geometries. 

VI. CONCLUSION AND FINAL REMARKS 

We have derived a Robinson-Trautman spacetime with 
minimally coupled free scalar field. We have shown that 
it has a singularity for all retarded times created by the 
divergence of the scalar field therein. This singularity is 
initially (with respect to retarded time) naked and only 
later becomes covered by the quasilocal horizon. Note 
that the energy momentum tensor of the free minimally 
coupled scalar field trivially satisfies null energy condition 
(as well as weak and strong ones) and the naked singular¬ 
ity at the beginning of the evolution is probably caused 
by a slow buildup of effective energy density caused by 
the scalar field at the singularity position which is enough 
to form the singularity but not enough to envelop it in 
horizon initially. This behaviour suggests similarity with 
the appearance of a naked curvature singularity in Vaidya 
spacetime with linear mass function. The naked singu¬ 
larity appears there initially depending on the speed of 
growth of mass [3l| and later becomes covered by hori¬ 


zon as well. From the properties of both null congru¬ 
ences orthogonal to the horizon we deduced that we are 
dealing with past horizon which is natural for standard 
(retarded) form of Robinson-Trautman spacetime. 

Our solution is asymptotically flat, contains a black 
hole (at least in the later stage of development) and has a 
scalar field, so one is naturally interested in its connection 
with the no-hair theorems (see [32j ] for current review). 
As recently shown [33|; for stationary black hole space- 
times there are no scalar hairs (even for time-dependent 
scalar field) which means that the dynamical nature of 
Robinson-Trautman family is truly needed for our solu¬ 
tion to be feasible. Also, we have shown that the scalar 
field vanishes outside the black hole in infinite retarded 
time limit when the geometry settles down to the final 
state — Schwarzschild black hole. 

Finally, we have proved that our geometry is of al¬ 
gebraic type II (the most general type for vacuum 
Robinson-Trautman spacetimes) and if we restrict to 
spherically symmetric case it is of type D. However, the 
type N subcase is not possible for our solution. 

Acknowledgments 

This work was supported by grant GACR 14-37086G. 


[1] I. Robinson and A. Trautman, Phys. Rev. Lett. 4, 431 
(1960). 

[2] I. Robinson and A. Trautman, Proc. Roy. Soc. Lond. 
A265 , 463 (1962). 

[3] H. Stephani, D. Kramer, M.A.H. MacCallum, C. Hoense- 
laers and E. Herlt, Exact Solutions of the Einstein’s Field 
Equations, 2nd edn (CUPress Cambridge, 2002). 

[4] P.T. Chrusciel, Commun. Math. Phys. 137, 289 (1991). 

[5] P.T. Chrusciel, Proc. Roy. Soc. Lond. A436, 299 (1992). 

[6] P.T. Chrusciel and D.B. Singleton, Commun. Math. 
Phys. 147, 137 (1992). 

[7] J. Bicak and J. Podolsky, Phys. Rev. D 52, 887 (1995). 

[8] J. Bicak and J. Podolsky, Phys. Rev. D 55, 1985 (1997). 

[9] J. Bicak, Z. Perjes, Class. Quantum Grav. 4, 595 (1987). 

[10] J. Podolsky and O. Svitek, Phys. Rev. D 71, 124001 

(2005). 

[11] J. Podolsky and M. Ortaggio, Class. Quant. Grav. 23, 
5785 (2006). 

[12] K.P. Tod, Class. Quantum Grav. 6, 1159 (1989). 

[13] E.W.M. Chow and A.W.C. Lun, J. Austr. Math. Soc. B 
41, 217 (1999). 

[14] J. Podolsky and O. Svitek, Phys. Rev. D 80, 124042 
(2009). 

[15] L. Rezzolla, R. P. Macedo and J. L. Jaramillo, Phys. Rev. 
Lett. 104, 221101 (2010). 

[16] O. Svitek, Phys. Rev. D 84 044027 (2011). 

[17] M. Ortaggio, J. Podolsky, M. Zofka, JHEP 1502 045 
(2015). 


[18] R. Giiven, E. Yorrik, Pliys. Rev. D 54 6413 (1996). 

[19] A. Anabalon, H. Maeda, Phys.Rev. D 81 041501 (2010). 

[20] J. B. Griffiths and J. Podolsky, Exact Space-Times 
in Einstein’s General Relativity, (Cambridge University 
Press, Cambridge, England, 2009). 

[21] S.W. Hawking and G.F.R. Ellis, The large scale structure 
of space-time (CUPress Cambridge, 1975). 

[22] S.A. Hayward, Phys. Rev. D 49, 6467 (1994). 

[23] A. Ashtekar and B. Krishnan, Phys. Rev. D 68, 104030 
(2003). 

[24] T-H. Kuo, Nonlinear Anal. 63, e427 (2005). 

[25] A. Besse, Einstein Manifolds (Springer Berlin, 1987). 

[26] L. Boccardo, F. Murat and J.P. Puel, Annali della Scuola 
Normale Superiore di Pisa (Classe di Scienze) Ser. 4 vol. 
11 no. 2, 213 (1984). 

[27] D. Gilbarg and N.S. Trudinger, Elliptic Partial Differen¬ 
tial Equations of Second Order (Springer Berlin, 2001). 

[28] E.g. supersolution satisfies Au + < F(x,u + ,Vu + ). 

[29] E. Zakhary, K. T. Vu, J. Carminati, Gen. Rel. Grav. 35, 
1223 (2003). 

[30] R. Penrose, W. Rindler, Spinors and Space-Time, Vol. 2 
(Cambridge University Press, 1990). 

[31] B. Waugh, K. Lake, Phys. Rev. D 34 2978 (1986). 

[32] C.A.R. Herdeiro, E. Radu [arXiv: 1504.08209]. 

]33] A.A.H. Graham, R. Jha, Phys. Rev. D 90 041501 
041501 (R). 



